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ABSTRACT 



Context. A search programme for pulsating subdwarf B stars was conducted with the Nordic Optical Telescope on La 
Palma over 59 nights between 1999 and 2009. 

Aims. The purpose of the programme was to significantly extend the number of rapidly pulsating sdB stars to better 
understand the properties of this new group of variable compact stars. 

Methods. Candidates were selected initially from the HS and HE surveys, but were supplemented with additional objects 
from other surveys. Short sequences of time-series photometry were made on the candidates to determine the presence 
of rapid pulsations. 

Results. In total twenty new pulsators were found in this survey, most of which have already been published and some 
extensively studied. We present four new short period pulsators, bringing the total of such pulsators up to 49. We also 
give limits on pulsation amplitudes for 285 objects with no obvious periodic variations,* summarise the results of the 
survey, and provide improved physical parameters on the composite pulsators for which only preliminary estimates 
were published earlier. 

Key words, subdwarfs - surveys - stars: oscillations - stars: individual: HE 2151-1001, HS 2125+1105, PC 1033+201, 
HE 1450-0957 



1. Introduction 

Hot subluminous stars are considered to be extreme hori- 
zontal branch (EHB) stars or closely related objects, with 
effective temperatures T e g ~ 20-35 kK. The EHB models 
imply that they are core helium burning objects with an ex- 
tremely th in (M env < 0.02M(^) inert hyd rogen dominated 



envelope (|Heberlll986t ISaffer et al.lll994D . This structure 
prevents them from ascending the asymptotic giant branch 
(AGB), and they must evolve instead towards higher tem- 
peratures and surface gravities after their core helium is 
exhausted. Thus, an sdB star crosses the hotter sdO do- 
main before reac hing degeneracy and cooling as a normal 
white dwarf star ()Dorman et al.lfl993[ ). However, important 
questions remain concerning the exact evolutionary paths 
and the appropriate timescales. 

How they evolve to the EHB configuration is contro- 
versial. The problem is how the mass loss mechanism in 



* The complete list with data on of 285 objects (Table [S| 
is only available in electronic form at the CDS via anony- 
mous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via 
|http : //cdsweb . u r s trasbg . f r/cgi-bin/qcat ? J/A+A/ ???/??? | 



the progenitor manages to remove all but a tiny fraction 
of the hydrogen envelope at about the same time as the 
core has attained the mass (~0.47M Q ) required for the 
He flash. About half of the sdB stars reside in close binary 
systems , with either a white dwarf or an M-dwarf as com- 
panion (|Maxted et aJj|2001t iNapiwotzki et al.ll2004f ). and a 
significant fraction of the rest are in wider binaries with 
a main sequence F-K star as companion. Therefore mass 
transfer in close binary evolution must be an important 
evolutionary pathway le ading to the formation of sdB stars 
(|Han et alj " l2002L I2003H. A recen t review of hot subdwarf 
stars is provided by iHeber] (|2009f ) . 

The discovery of multimode pulsations in sdBs opened 
an attractive opportunity of probing their interiors with 
seismological methods. The properties of the sdB pul- 
sator s (sdBVs or V36 1 Hy J] stars after the proto- 
type; iKilkennv et al.lll997[ ) are characterised by relatively 
short pulsation periods (~1 to lOmin) and low pul- 
sation amplitudes attributed to low order pressure (p) 



1 Translations between the old survey names and the variable 
star names used here are provided in Table [9] 
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modes. Most V361Hya stars have been found with pul- 
sation amplitudes of around ten millimodulation am- 
plitudes (mmg0), although a few objects show a main 
peak amplitude up to ~60 mma (V338 Ser. iKilkennv et ahl 
119991: Balloon 090100001, lOreiro et al.l |200j), while the 
lowest level pulsato r discovered up to now is LMDra 
(jSilvotti et al.l [20001 ) . with no peaks above 2mma. The 
number of detected periods varies from a single one to more 
than fifty, strongly correlated with the accuracy of the mea- 
surements. 

In 2001, some sdBs were disco vered to show long -period 
(~lh) photometric modulations (|Green et al.ll2003[ ). which 
were interpreted as high radial order gravity (g) modes. 
These stars are referred to as long period sdBV stars or 
V1093Her stars after the prototype. V1093Her stars are 
cooler than their short-period counterparts, although their 
instability regions appear to overlap around T c g ^29 kK, 
where they can display simultaneous short and long period 
pulsations. These stars are referre d to as hybrid sdB Vs or 
DWLyn stars after the prototype (ISchuh et al. 2006f). 

Non-adiabatic computations bv lCharpinet et al.l ( 19971 ) 
predicted excited low degree (£), low radial order (n) p- 
modes in sdB models, driven by an iron opacity bump, in 
the temperature range 29-37kK. This mechanism is in- 
efficient at solar metallicity, but gravitational settling and 
radiative levitation can work together to locally enhance 
metals in a driving zone in the envelope. This k mechanism 
has been successfully invoked to explain both the p-mode 
pulsations in V361 Hy a stars and the q- m ode pulsations 
in V1 093Her stars (jFontaine et al.ll2003D . Uefferv fc Saiol 
(|2006h were able to expand the instability regions by includ- 
ing opacities for nickel, to the point where the p-mode and 
g-mode domains overlap, thereby explaining the presence 
of hybrid oscillations. However, it is still not understood 
why most stars in this temperature range do not appear to 
vary. 

As the number of known short-period and hybrid sdBV 
stars have increased, several peculiar features are emerging: 

— Mode density: The number of pulsation frequencies de- 
tected vary from a sing le one to more than 50. LS Dra 
(|0stensen et al.l[200TarT was confirmed to be monoperi- 
odic well below the millimagnitude lev el from 127 hours 
of observations bv iReed et al.1 (|2007bl ). For others, de- 
tailed follow-up observations have increased the number 
of detected frequencies to more than 50. 

— Period grouping: Pulsation periods for the whole 
V361Hya group lie between 90 and 800 s, but distinc- 
tions appear to be present between the stars in the 
shorter and longer period groups. All the sdBVs with 
periods between 300 and 400s except one are DWLyn 
type hybrid pulsators. The odd one out is KLUMa, 
which also stands out as th e only binary of the p ulsators 
with periods in this range (|Q'Toole et al.ll2004l) . 

— Amplitude variations: Amplitudes are found to change 
with time in many cases. However, it is unclear whether 
this is due to true a mplitude variati ons or beating of 
closely spaced modes (|Kilkenny||2010l ). 

— Dominant modes: According to their amplitudes the 
stars can be grouped into high amplitude pulsators, for 



which a mode was observed to exceed 30 mmag at least 
once, and low-amplitude ones for which such dominant 
modes were never reported. All the high amplitude pul- 
sators are DWLyn type, except the unique V338Ser, 
which lies well above the canonical EHB region. 

The V361 Hya stars with very few exited modes are not 
suitable for the period-matching technique in asteroseismol- 
ogy, since too many models can be found that adequately 
match their sparse periodograms. However, if a mode has 
sufficient amplitude, the degree t can be determined from 
spectroscopy, substantially constraining the asteroseismic 
solution. For low-amplitude V361 Hya stars with sufficiently 
rich spectra, however, asteroseismology has met with great 
success, as the most fundamental parameter of a star - its 
mass - has been derived for ten such pulsators a long with 
their fractional hydrog en envelope masses ( see 0stenscn 
120091 for a review, and iRandall et aLll2009bl for the most 
recent result). 

Nevertheless, sdBVs with sparse pulsation spectra are 
important as they may allow the measurement of light 
travel time variations or period changes due to stellar evo- 
lution, if the modes are sufficien tly stable, as exemp lified 
by the planet-hosting V391Peg (ISilvqtti et al.1l2007tl . dis- 



2 The mma units are 10 -3 of the Fourier amplitudes of a light- 
curve in normalised intensity units. One mma translates to a 
peak-to-peak amplitude of two millimodulation intensity (mmi) 
units in the light curve. 



covered by our survey ( 0stensc n et al.ll2001bD . Hence, pho- 
tometric surveys to discover new V361 Hya stars are still 
rewarding. 

The theoretical location of the sdBV instability strip 
in the T e ff-log<7 plane is quite well established. However, 
no more than one out of ten sdB stars located within 
the instability region are actually found to pulsate. This 
situation is quite different from that of white dwarf pul- 
sators of the ZZ Ceti type, for which evidence is growing 
that all white dwarfs in their instability strip are pulsating 
([Gianninas et al.l 12007). For sdBVs it has been suggested 
that younger EHB stars may not have enough iron accu- 
mulated in the driving region for the K-mechanism to be 
efficient enough to dri ve pulsations, but tim e-dependent 
diffusion calculations bv lFontaine et all ([2006) have demon- 
strated that sufficient iron for pulsations to occur accumu- 
lates after a few hundred thousand years, and after 1 Myr 
no further accumulation is achieved. Compared to the EHB 
lifetime of 100-150 Myr this is not significant. An interest- 
ing speculation was made bv lJefferv fc Saiol (j2007l) , that the 
iron group element enhancements may be disrupted by the 
atmospheric motions as pulsations build up to some level. 
They note that since p-modes mostly involve vertical mo- 
tion, while g-modes are dominated by horizontal motion, it 
is possible that p-modes are more effective at redistribut- 
ing the iron group elements out of the driving zone. This 
could explain the observation that most cool EHB stars 
are g-mode pu lsators while at th e hotter end of the branch 
most are not (|Green et al.l l2003h . It may also be that the 
shortage of pulsators is due to observational biases, either 
that the pulsations have too low amplitudes, or that they 
have too high £ to be easily detectable photometrically, an 
issue that needs further investigation. To this end it also 
is important to publish observational constraints on stars 
that were not found to be variable in order to guide future 
more sensitive surveys. 

Here we present the final results of a long-term pro- 
gramme to detect new short-period sdBVs using the Nordic 
Optical Telescope (not) on La Palma. 
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Table 1. Observation runs at the Nordic Optical Telescope 
contributing to this survey. 



ID 


Dates 


AT 


AT 


Observers 


Nl 


1999 July 19-23 


4.0 


13 


JES, RS, JMGP 


N2 


1999 Oct. 14-19 


5.0 


31 


JES, RH0 


N3 


2000 July 7-11 


4.0 


25 


RH0, RS 


N4 


2000 Oct. 4-8 


4.0 


38 


JES, RH0 


N5 


2001 July 25-29 


4.0 


22 


RHO, RS 


N6 


2002 Oct. 10-15 


0.5 


4 


JES a 


N7 


2004 June 4-8 


4.0 


44 


JES, RHO 


No 


2005 Feb. 15-20 


5.0 





JES b 


N8 


2006 Dec. 10-15 


0.3 


4 


JES a 


N9 


2007 Dec. 14-16 


3.0 


14 


RHO, JHT C 


XI 


2004 Oct. 15-17 


0.3 


4 


JES, RH0 d 


SI 


2004 May 21-24 


4.0 


20 


RO 


S2 


2005 May 30-June 1 


4.0 


15 


RO 


C Q 


nnn? f~\s*4- OA O A 

2UU/ Uct. zU— z4 


q n 
O.U 


1 1 


tiU 


^ A 

OA 


9fin& C!or\ 99 9^ 
ZUUo Dep. zz-zo 


o.u 


zu 


1 \ w 


S5 


2008 Tulv 26-30 

— ' ' V ' v..' • } Illy _■ \ ) • t V 7 


2.0 


67 


RH0 c 


S6 


2009 Jan. 24-27 


3.0 


5 


RO 


Tl 


2001 Feb. 24 


0.6 


7 


RH0 f 


T2 


2001 April 10, 13 


2.0 


24 


RH0 e 


T3 


2002 Sep. 13-14 


0.2 


4 


RHO, JHT f 


T4 


2003 Feb. 10-19 


0.5 


5 


RHO, JHT f 


T5 


2004 Jan. 13 


0.5 


9 


RHO, JHT f 


T6 


2005 May 5 


0.7 


16 


RH0 f 


T7 


2006 May 9 


0.7 


15 


RH0 f 


T8 


2008 Feb 24, 26 


0.2 


2 


RH0 f 


ST 


2006 June 10 


0.5 


2 


JMGP g 


Total 




59.0 


421 





The ID column gives our run identifier, used to refer to a partic- 
ular observation run in Table [5] and when discussing particular 
observations. N n gives the number of nights used for this sur- 
vey programme. For nights allocated to the survey, the number 
includes all nights, regardless of whether any observations were 
actually made or not. For nights where this programme served 
as a backup for another, only the fraction of nights actually used 
for this programme is listed. The Nt column shows the number 
of targets observed. Note that reobservations and observations 
of sdO stars (to be analysed in a separate paper) are included 
in these totals. The notes indicate 
a Run primarily dedicated to an unrelated target. 
b Run lost due to snow and ice. 
c Run primarily for follow-up of our pulsators. 
d Xcov 24 multi-site campaign. 
° Run primarily dedicated to search for sdO stars. 
f Run primarily for technical investigations. 
g Spanish service time. 



2. The survey 

Shortly after the discove ry of the first pulsat i ng sdBs by 
the South African group (Kilkenny et al.l[l997t iKoen et al.l 
Il997t IStobie et alj|l997at lO'Donoghue et al.lll997l) . it was 
decided to initiate a programme to search for such pulsators 
using the 2.5 m NOT on La Palma. The principal aim was 
to extend the number of sdB pulsators, determine the frac- 
tion of variables within this class of stars, and explore the 
boundaries of their instability region. The South African 
group had at the time discovered and confirmed pulsations 
in 14 objects from a sample of around 600 spectroscopically 
confirmed hot subdwarf stars, and a Canadian team found 
four pulsators in 74 stars after selecting candidates based on 
T e ff (jBilleres et alJ[2002h . Since those samples were overlap- 



ping considerably, two of the pulsators were found in both 
samples. To avoid searching an already depleted sample, we 
drew our targets primarily from a new spectroscopic study 
targeting sdB star candidates fr om the Hamburg-Schmidt 
(HS) survey ()Hagen et al.l Il995h . although it was comple- 
mented with targets from other surveys as explained below. 

The complete list of observational runs fully or partly 
used for this survey is presented in Table [TJ In Table [5] 
we provide coordinates and observational limits on all sdB 
stars checked for variability during the survey, except those 
that were found to pulsate or vary due to eclipses and re- 
flection effects. The eclipsing and reflection variables are 
published in separate papers, and the 20 pulsators are listed 
in Table 13 together with pulsators found by other surveys. 

2.1. Chronology and results published to date 

The rapid pulsations in sdBVs are a challenge to observe 
with standard CCD instrumentation due to the pulsation 
periods being comparable to the typical readout times of a 
full CCD frame. An essential part of this project was the 
design of a system capable of efficien t windowed CCD pho- 
tometry (|0stensen fc Solheiml [20001. A dedicated camera 
was built for this purpose, the Troms0 CCD Photometer 
(TCP), and the control system implemented for this cam- 
era was adopted to both the hirac and ALFOSC cameras at 
the NOT. 

The first search run was made in the summer of 1999 
(c.f. Table H]), using the Troms0- Texas Photometer (TTP) 
and revealing the first sdB pulsator from this programme 
(jSilvotti et all l2000l Paper i). Already the same autumn 
the CCD system was operational, and the improved detec- 
tion efficiency allowed us to reach acceptable noise levels 
in a much shorter time, thereby resulting in the discov- 
ery of three new pulsators in only five nights of observa- 
tions (|0stensen et al.l l2001bl Paper n) . Two observation 
runs in 2000 led to th e discovery of four more pulsators 
()0stensen et aljfeOOlal Paper in). At this time we started 
to run out of targets from the HS survey at some RA's, 
and supplemented our list with stars from the PG survey 
(iGreen et al.l fl986h This was a concern, as many of the 
known sdB stars had already been targeted by the South 
African group, and we expected our detection efficiency to 
drop since we avoided targeting stars already reported to 
pulsate. Still, the sample was not completely depleted, and 
during the summer 2001 run we discovered pulsations in 
QQVir (PG 1325+101) and EPPsc which originated from 
the HS sample (HS 2303+0 152) but also appear in the PG 
catalogue as PG 2303+019 (jSilvotti et al.ll2002bl Paper iv). 

To get a sample that was both sufficiently large to prop- 
erly determine the extent of the sdB instability region and 
the fraction of pulsators across this region, and at the same 
time undepleted, it was decided to sup plement our tar- 
get li sts with stars from the SDSS survey (jStoughton et al.l 
2002). The first three SDSS stars were observed in October 
2002, and one of them, SDSS J171722.1 0+580559.9 or 
J171 7+5805 for short, was found to pulsate (Solh eim et al.l 
12004 Paper v). 

In spite of the initial success of our programme, the 
Nordic time-allocation committee started turning down 
our proposals for further search time, and the last regu- 
lar search run on the original programme was allocated 
for February 2005, but was completely lost due to bad 
weather. The last successful search run was made in 
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June 2004, resulting in three new pulsators. These are 
PG 1419+081, SDSS J144514.93+000249.0 (J1445+0002) 
and SDSS J164214. 21+425234.0 (J1642+4552). Only brief 
sequences w ere made on these ob j e cts, a nd the results were 
presented in lSolheim fc 0stensenl (|2006l Paper vi). 

At this time a collaborative search was initiated with 
a Spanish team (RO, AU, FPH, CRL), which permitted 
additional access to the NOT through Spanish time. With 
the Spanish group, we were able to explore a la rger num- 
ber of stars drawn from the Subdwarf Database ()0stensenl 
12004 12006|) . The first pulsator discovered with the NOT on 
Spanish time was another star drawn from the SDSS sample 
PG 1657 +416 (SDSS J 1 65841 .83+413115.6), and was pre- 
sented by lOreiro et al.l (|2007l . Paper vn) . Observation runs 
allocated on Nordic time are indicated by run numbers pre- 
fixed by N, and on Spanish time by S in Tabic [1] 

Nordic time still continued to be granted for follow-up 
programmes on various pulsating stars, including one of 
the most interesting pulsators discovered during this pro- 
gramme, V391 Peg, which was demonstrated to be host to a 
planet bv lSilvotti et al.l (|2007l ). by using the orbital period 
modulation introduced on the two main pulsation modes. 
A few unexplored pulsator candidates were observed dur- 
ing these runs (N6, N8, N9 in Table [T]), and also during the 
WET Xcov 24 campa ign on the pulsating sdB star UY Sex 
([Vuckovic et al.ll2006l run XI), whenever the prime target 
was not reachable. Upgrades and tests on the fast CCD 
photometry system also gave several opportunities to ob- 
serve some additional targets during technical time at the 
not (listed as Tl to T8 in Table [T]). One new pulsator, 
PG 1033+201, was found during run T7. 

The limitations on telescope time from 2004 onward 
forced us to abandon the faint SDSS targets. Instead, we 
focused on detecting more pulsators in the cool region 
where the particularly interesting objects V391Peg and 
Balloon 090100001 had been found. The discove ry of hy- 
brid pulsations in DW Lyn by ISchuh et al.l ([2006) provided 
further encouragement for increasing the sample in this re- 
gion of the instability strip. Although no more such cool 
pulsators were detected, a total of five more pulsators in 
the hotter end of the instability strip were found. One, 
from the BG sample (Bok-Green iGreen et al.l [2008'). was 
recently published in lOreiro et al.l (|2009l . Paper vm), and 
four are presented for the first time in this paper. 

In addition to these 20 pulsators our photometric ob- 
servations have revealed several of the relatively rare sdB 
stars with M-dwarf companions from their strong reflec- 
tion effect. The discovery of an eclipsing sdB+dM sys- 
tem, HS 07 05+6700, first detected during run N4 was pre- 
sented by iDrechsel et al.l (|200lh . and the non-eclipsing 
sdB+dM system HS 2333+39 2 7 disc overed during run N2 
was presented bv lHeber et al.l ([2004). Two more sdBs with 
dM companion s were found during run S2: the eclipsing 
HS 2231+2441 (lOstensen et al.ll2007h . and HS 2043+0615 
()0stensen et all hi prep J ). 

Our survey was not originally designed to detect long 
period variability, since we aimed for short photometric se- 
quences with high S/N. However, the large photometric 
amplitudes of sdB+dM binaries produce a strong trend in 
the differential photometry, of which we saw several exam- 
ples. However, since follow-up of such long-period variable 
objects can easily be done with smaller telescopes using 
long integration times, no follow-up was attempted on these 
stars with the NOT. The long period g-mode sdBVs also pul- 



Table 2. Analysis for eight stars from the HS survey, not 
included in lEdelmann et all ([20031 ) . 



Name 


m pg 


Class 


T c ff 


logs 


logy 




mag 




kK 


dex 




Observed with CA/twin 










HS 1733+4540 


15.1 


sdOB 


37.5 


5.6 


-3.0 


HS 1909+7004 


15.4 


sdB 


36.5 


5.4 


-2.9 


Spectra courtesy 


of AEuilar-Sanchez (1998) 




HS 0026+0439 1 


14.8 


sdB 


36.7 


5.8 


-3.0 


HS 0042+0927 1 


16.0 


sdB 


30.2 


5.6 


-2.9 


HS 2254+0640 1 


15.2 


sdB 


29.1 


5.6 


-2.7 


HS 2320+0840 1 ' 2 


14.6 


sdB+ 


30.0 


5.8 


-2.6 


HS 2323+0459 1 


14.8 


sdB 


29.3 


5.4 


-2.8 


HS 2334+0144 


15.5 


sdB 


30.8 


5.6 


-2.6 



m pg gives the photographic magnitude from the HS survey 
plates, 'Class' our spectroscopic classification, followed by the 
effective temperature (T e ff), surface gravity (logg), and at- 
mospheric helium abundance (log y = log (Aho/Ah)), from our 
model fits. The sdB+ classification indicates a composite spec- 
trum. The notes on the names i mply: 

1 Also found in the PG survey (IGreen et alJl!986f). 

2 Also found in the KUV survey ijNoguchi et al.lll980f ). 



sate on time-scales around one hour, but with much lower 
amplitudes. As they had yet to be discovered at the start of 
our survey, our observation strategy was not designed with 
these objects in mind. 

2.2. The HS/ HE/SPY samples 

When our programme was initiated, the main source of 
targets was the HS survey, but we also included targets 
fro m the equatorial part of the Hamburg-ESO (HE) sur- 
vey (jWisotzki et al.l[l996T ). Follow-up spe ctroscopy made at 
Calar Alto with the tw in spectrograph (jHeber et al.lll999l: 
lEdelmann et all 120031 ) allowed us to preselect candidates 
with effective temperatures and gravitie s in the domain 
predi cted for the pulsational instability (jCharpinet et al.l 
119971 ). This provided us with a sample of pulsator can- 
didates that was both tuned to the instability region 
and relatively undepleted by previous surveys. The lat- 
ter was particularly important since the South African 
group had already explored a considerable fraction of 
the sdB stars known in the literature (and compiled in 
the C atalogue of spectrosco pically identified hot subdwarf 
stars, iKilkennv et all [l988l However, it has turned out 
that the HE survey has a significant overlap with unpub- 
lishe d parts of the Edi nburgh-Cape (EC) blue object sur- 
vey (jStobie et alJll997bh , which was one of the main sources 
surveyed by the South African group. This can explain why 
out of 80 stars from the HS survey and 20 stars from the 
H E survey observed during the initial five runs (as reported 
in lPaperVf l. pulsations were found in eight HS stars and in 
none of the HE stars. Due to the difficulty of establishing 
accurate effective temperatures and gravities for sdB+F- 
K stars, these were not included in our initial sample, and 
were only added later when the priority sample started to 
become depleted. Thus, unlike the South African group, 
who found that most of their first pulsators were in binary 
systems, due to the selection biases involved only one of 
our initial 10 were. Note that another composite, our first 
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pulsator LM Dra, is in a well-separated visual binary with 
an F3 companion, and therefore not considered to be an 
sdB+F system since the stars cannot have interacted dur- 
ing their evolution. 

Only eight stars from the HS sample observed during 
this survey have not yet been published elsewhere. Their 
coordinates are listed in the on-line table together with 
those of all the other targets (Table [SJ, and their phys- 
ical parameters are listed in Table [5J We computed the 
stellar atmospheric parameters (effective temperature, sur- 
face gravity, and photospheric helium abundance) by fitting 
model atmosphere grids to the hydrogen and helium lines 
visible in the spectra. The procedure us e d for this fitting 
was the same as that of lEdelmann et al.l (l2003h but using 
only the LTE models described in iHeber et al.l {2000), as 
the purpose of our modelling was only to establish which 
stars are roughly within the sdB instability region rather 
than to establish precise physical parameters. 

The parameters of HS 1909+7004 are uncertain since 
the spectrum is of low resolution and the fit is quite 
poor, but the fits to the other stars appear to be quite 
good. HS 2320+0840 has a very strong companion that con- 
tributes clearly visible line features in the spectrum. For 
this reason, we performed a spectral decompositio n sub - 
tracting a main sequence model from iMunari et~aTl (l2005h 
with T e fj = 6.5 kK, log g — 4 and solar metallicity, which we 
found to contribute about 40% of the flux at 5 000 A. 

A total of 46 objects from the HE survey were observed. 
Fo r 23 of the s e the physical parameters were determined 
by lEdelmannl (|2003|) , and none of these were found to pul- 
sate. Cross-referencing them with the literature reveals that 
only nine are new, the rest are subdwarfs with identifica- 
tions from the Montreal- Cambridge- Tololo (MCT) Survey 
(lLamontagne et ah! l2000t) . or earlier surveys compiled by 
Kilkenny et al.l (|1988l ). 

A number of white dwarf candidates fr om the HE survey 
were observed as part o f the SPY survey (iNapiwotzki et all 
120011) and published bv lLisker et"aL1 (|2005TT where physical 
parameters are provided also for the composite objects. 23 
HE stars from the SPY sample were observed, and two were 
found to pulsate. Also, four HS stars from this sample were 
observed, and one more pulsator was found. Thes e three 
new variables are presented in Sect. [51 iLisker et al.l ([2005) 
also present a list of sdB stars formerly misclassified as 
white dwarfs. We observed four of these, but did not detect 
signs of variability^ 

2.3. The SDSS sample 

More than 16000 spectra of UV excess stars were down- 
loaded from the SDSS survey (jStoughton et al.i r2002) up to 
data release 6 (DR6), and the hot subdwarfs were iden- 
tified by inspection. The intermediate resolution (~1 A) 
SDSS spectra cover all wavelengths from 3 800 to 9 000 A, 
ideal for classification purposes. We identified about 900 
objects as hot subdwarf stars and classified them into sdB, 
sdO, sdOB, He-sdB, He-sdO, and He-sdOB, depending on 
whether hydrogen or helium dominates the spectrum, and 
if they contain lines from Hei, Hen or both. 318 of the sdB 
stars were classified as single sdB stars and 165 as sdB+F-K 



Table 3. The stars from the SDSS survey observed during 
this programme, excluding the pulsators listed in Table [4] 



Name 


g' 


Class 




logs 


logy 




mag 




kK 


dex 


dex 


PB5916 1 


15.3 


sdB 


30.1 


5.7 


-2.3 


PG0812+482 2 


15.0 


sdB 


25.0 


5.4 


-3.0 


PG 0826+480 2 


15.8 


sdB 


26.1 


5.5 


-3.1 


PG 1100-008 2 


16.3 


sdB+ 


29.2 


5.5 


-2.6+ 


PG 1136-003 1 ' 2 


14.2 


sdB 


31.6 


5.6 


-2.7 


PG 1249+028 1 ' 2 ' 3 


15.4 


sdB 


30.6 


5.6 


-3.0 


J12596-0039 4 


16.8 


sdOB+ 


34.4 


5.6 


-2.2+ 


PG 1315+013 1 ' 2 


16.6 


sdB 


26.6 


5.3 


-2.3 


J13320+6733 4 


17.0 


sdOB 


36.1 


6.0 


-1.1 


J13516+0234 1 


17.1 


sdB+ 


34.7 


5.3 


-3.0+ 


PG 1403+019 1 ' 2 


15.6 


sdB 


27.7 


5.5 


-1.9 


J14086+6531 1 


17.4 


sdB 


30.4 


5.6 


-3.0 


J14236+0149 1 


17.1 


sdB 


29.5 


5.7 


-2.3 


PG 1422+035 2 


16.2 


sdOB 


35.1 


5.8 


-1.5 


J15107+0409 4 


16.9 


sdOB 


34.7 


5.7 


-1.6 


J15131+0114 1 


17.1 


sdB 


28.0 


5.6 


-3.0 


EGGR 491 1 ' 5 


16.8 


sdB 


30.6 


5.6 


-3.0 


SBSS 1544+568 1 


16.8 


sdB 


27.2 


5.2 


-2.0 


J15481-0049 1 


16.4 


sdB+ 


32.7 


5.6 


-3.0+ 


J15513+0649 


15.9 


sdOB+ 


34.3 


5.0 


-2.1+ 


J15564+0113 1 


16.0 


sdB 


30.7 


5.7 


-3.1 


J16026-0012 1 


17.5 


sdB 


34.0 


5.8 


-2.1 


J16165-0038 1 


16.7 


sdOB 


36.0 


5.7 


-1.8 


J16331+0032 1 


16.9 


sdB 


31.7 


5.8 


-3.3 


J16347-0053 1 


16.9 


sdB 


29.1 


5.5 


-2.7 


J16420+4403 1 


16.7 


sdB 


28.8 


5.2 


-2.8 


J16443+4523 1 


17.1 


sdB 


33.1 


5.8 


-2.0 


PG 1653+633 1 ' 2 


15.9 


sdOB 


36.1 


5.9 


-1.5 


J17144+6147 1 


16.6 


sdB+ 


33.4 


5.6 


-3.4+ 


SBSS1715+556 1 ' 6 


16.9 


sdB 


30.3 


5.4 


-3.2 


J17165+5751 4 


18.0 


sdOB 


34.9 


5.8 


-0.8 


J20440-0511 1 


17.3 


sdB 


29.9 


5.4 


-3.2 


J20457-0543 1 


17.8 


sdB 


34.9 


5.4 


-1.7 


J20573-0626 1 


16.5 


sdB 


28.0 


5.1 


-2.5 


J21531-0719 1 


16.9 


sdB 


32.0 


5.9 


-2.0 



The columns are as in Tabled except that Sloan g' magnitudes 
are provided provided as indicator of the signal level. A plus 
sign after the class marks the presence of a cool companion from 
features in the spectrum, the slope of the spectrum or z-band 

excess. The notes indicate: 

Correctly identified as sdB stars in Eiscnst ein et all lf2 006h 
Classified as hot subdwarf stars in lGreen et all (|1986f l. 
Object misnamed in t he PG catalog; should be PG 1249-027. 

Listed as sdO stars in Eisenst ein et all (|2006[). 

Misclassified as a white dw arf in McCook fc Sionl lll987f) . 
Identified as an sdB star in lSteparnari*^t^!Tl|2OTll h 
Parameters for composites are less certain than for singles. 



3 Actually, one more of the fourteen stars in this list of mis- 
classified WDs is included in our sample: TonS 155 appears in 



lEdelmannl [|2003D as HE 0021-2326. 



binaries (34%) based on spectroscopic signatures and their 
position in (u' — g',g' — r') colour-colour space. A similar 
fraction was found for the sdOB stars; 78 single and 35 bi- 
naries (31%). All these stars were fitted to model spectra 
to determine their position in the T c g /log g plane, using 
the Balmer lines Ha to Hij, as well as the most promi- 
nent helium lines at 4472 A and 4026 A for the sdB stars, 
and including Hell 4686 for the sdOBs. The formal fit- 
ting errors are about 250 K for T c s, 0.05 dex for \ogg, and 
0.1 dex for log y = log (A^ho/A^h) at a g' magnitude around 
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Table 4. Spectral decomposition of the binary pulsators from the SDSS sample. 



Companion signature Primary Secondary 
Name Can CH Mgi Nai E(B — V) / 50 oo T eS log ff logy / 50 oo T cff logg [M/H] 
H+K G b D mag kK dex kK dex 



J1717+5805 


v 7 




V 


V 


0.025 


PG 1657+416 




V 


V 




0.021 


PG 1419+081 


V 






si 


0.028 


J 1445+0002 










0.11 


J1642+4552 


V 


V 


V 


V 


0.011 



16.0, but increase by a factor of two for the fainter stars 
between g' = 17.5 and 18. For the stars with F-G compan- 
ions the uncertainties are larger, but not so large that the 
parameters are not useful for establishing whether or not 
the stars are located within the instability region. We note 
also that for the hottest st ars in the sample, NLTE effects 
start to become significant (Napiwotzki 1997), and that the 
LTE models can underestimate the temperatures by up to 
1000 K in some cases. 

Only 40 stars from this large sample were observed due 
to observing time restrictions. Five were found to pulsate 
and are listed in Table S] with additional details in Table GO 
The remaining 35 objects are listed in Table [3l with the 
physical parameters we derived from the SDSS archive spec- 
tra. The sample includes eleven spectroscopic binaries, and 
curiously all five pulsators found from the SDSS sample 
were among these. The physical parameters for the bina- 
ries have considerable systematic shifts depending on the 
contribution from the main sequence companion. For this 
reason we have made a spectral decomposition of the five 
pulsators in order to place them reliably in the T e ff/log<? 
plane. The details of the spectroscopic decomposition are 
given in Table SJ The six non-pulsating sdB+F-K systems 
were fitted as if they were single stars. 

It is non-trivial to estimate the spectroscopic parame- 
ters of the companion, but a reasonable compromise can 
usually be made that approximates the spectral contribu- 
tions from the main sequence companion to the Can H 
and K lines as well as the g-band and the Mg I lines. For the 
spectral decomposition we used the SDSS spectra processed 
with the DR7 pipeline, which have a significantly improved 
flux calibration compared with earlier data releases. After 
dereddening the spectru m by the E(B — V) re ddening co- 
efficients provided by the ISchleeel et al.l (|l998h dust maps, 
the flux was fitted with one of our subdwarf model spectra 
and a main sequence model spectrum from iMunari et al.l 
(|2005l) . The parameters of the subdwarf are then refitted 
to the observed spectrum after subtraction of the main se- 
quence model, and the procedure was iterated until it con- 
verged. Note that the reddening coefficients are for lines of 
sight to infinity and that we have no guarantee that the 
particular subdwarf under consideration is not in front of 
the bulk of the dust. There is also a significant uncertainty 
in the E(B — V) valu es due to the limited resolution of the 
ISchleeel et all (fl998l) maps and the lumpiness of the inter- 
stellar medium. But due to the low resolution and limited 
signal in the SDSS spectra, it is hard to make a complete 
disentanglement of the spectra while leaving the reddening 
as a free parameter, so we will stick with the dust map val- 
ues. In the case of J1445+0002, which has the highest red- 
dening coefficient of the five, it is quite clear that the Call 
lines are too narrow to originate from the main sequence 



0.90 


34.4 


5.75 


-1.8 


0.10 


5.0 


4.0 


+0.0 


0.74 


32.2 


5.80 


-2.0 


0.25 


5.5 


4.0 


+0.0 


0.86 


33.3 


5.85 


-1.8 


0.14 


5.5 


4.0 


+0.0 


0.83 


35.9 


5.75 


-1.6 


0.17 


6.0 


4.0 


+0.0 


0.53 


32.4 


5.80 


-2.0 


0.47 


6.0 


4.0 


-1.5 



companion, and no solution can be found that removes the 
Call contribution from the subdwarf spectrum. Therefore, 
the Ca II signature is mostly interstellar in origin, consistent 
with the high reddening value. The strongest companion by 
far is that of J1642+4552, in which it contributes almost 
half the light at A = 5000 A (/ 500 o in Table SJ). The system 
is also in a region of very low extinction, reducing the ambi- 
guity of the spectral decomposition. The strong companion 
and good signal in this particular SDSS spectrum allow us 
to make some further constraints on the companion. It is 
clear that a solar metallicity makes a far stronger contri- 
bution to the metal bands than allowed by the observed 
spectrum. A better fit is achieved with a spectrum metal 
depleted relatively to solar by [M/H] = —1.5. 

The temperatures derived here are somewhat lower than 
the temperatures we originally used to select the stars, but 
not enough to make a diffe rence with respect to the sample 
selection, as we suspected (jPaper VI J ) . For J1445+0002 the 
temperature drops from 37.6 to a more reasonable 35.9 K. 
It is suspicious that the derived log g values for all the five 
composite stars are so similar (within ±0.1 dex), and might 
be due to degeneracies when fitting so many free parame- 
ters to spectra with rather low resolution. It has also been 
noted that more reliable estimates of the surface gravity are 
obtained when the spectra cover all the high order Balmer 
lines, a requirement the SDSS spectra do not satisfy. 

2.4. The Bok-Green sample 

Eleven stars in our survey were selected from a spectro- 
scopic campaign to obtain a large unbiased s ample of hot 
stars based on 2MASS colours, undertaken by iGreen et al.l 
(|2008l ) with the University of Arizona 2.3 m Bok telescope. 
Physical parameters were de rived from a su bsample of 89 
stars from this campaign bv IWinterl (|2006l Appendix C). 
None of the stars we observed have reliable spectroscopic 
classifications in the literature, but three have been listed as 
faint blue stars by earlier surveys (Ton 930, FBS 1133+754, 
and FBS 1224+780). However, none of these were known 
to be hot subdwarf stars, so we consider all eleven as 
belonging to an undepleted sample. One of the sdBs in 
the BG sample was found to pulsate, 2 M0415+01 54 = 
2MASS J04155016+0154209, as reported in lPaper vnj The 
remaining seven stars are listed with 2M designation in 
Table [3 

2.5. Literature 

To supplement our sample with targets to fill all RAs we 
added stars with temperatures and gravities published in 
the literature. Most of these additional objects are well- 
known subdwarfs from the PG survey and are included 
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in many studies. The largest sample used was that of 
ISaffer et al.l JT994) , where 24 targets were surveyed and one 
found to pulsate (LMDra). 

Seventeen PG st ars from the radial velocity survey of 
(120011) . and four stars from the HST study of 
were also included in our sample, but 



Maxted et al 



Heber et al 



none were found to be variable. 

Rough temperature estimates can also be made on the 
basis of photometric measurements alon e, and severa l works 
provide such estimates. The survey of iBeers et all (|1992h 
includes temperatures computed from UBV colours, and 
we observed eight of these stars. A handful of c ompos - 
ite sdB+MS stars from the list of lAllard et all (j 1994b . 
where their estimate of the temperature for the primary lie 
within the instability region, were also observed. Another 
four star s were targe t ed ba sed on the temperature esti- 
mates of iBixler et al j ( [ 1991D. Simil a rly, w e observed stars 
from the sample of Moehler et al.l (|1990l ). after estimat- 
ing the temperatures based on their Stromgren photom- 
etry and the spectral features visible in their spectroscopic 
atlas. In total, about fifty stars in the sample were ob- 
served based on such photometric temperature estimates, 
and only one, PG 1419+081, was found to vary. The SDSS 
spectroscopy of this object was released after our discovery, 
allowing us to refine the original colour temperature esti- 
mate and correctly identify the star as an sdB+G5 compos- 
ite. We found tha t the colour temperature (T e g = 33.8 kK; 
IBeers et al.l [l992) is encouragingly close to our model fit 
value (T eff = 33.3kK). 

If we look at the complete sample from IBeers et al.l 
(|l992tl and consider only the stars classified as sdB with 
temperatures above 27.5 kK we are left with 23 stars, of 
which we have observed nine and found one variable. But 
EC 20338-1925 is also in this sample and was in our origi- 
nal list, although excluded after pulsations were discovered 
by th e South African team (published by Kilkenny et al.l 
2006, but preliminary results from the 1998 discovery were 
distributed much earlier) . A nother star in the sample is 
BPS CS 22169-1, analysed bv lGeier et all (I2009T) and found 
to have T P a = 39.1 kK, much hotter than the 33.9 kK es- 
timate by IBeers et al.1 (|1992l ). and hotter than the pre- 
dicted sdBV instability region. The star BPS CS 22890-94 
= PG 1525+024 was not observed s ince the colour tempera- 
ture estimate by IBeers et al.1 (|1992l ) placed it below our cut 
at T e ff = 22.9kK. However, recent SDSS spectroscopy gives 
us T e ff = 27.9 kK, which would have placed the object inside 
our cuts, further emphasising the uncertainties associated 
with using colours as a temperature estimator. 



3. Observations 

During our observations we systematically aimed at identi- 
fying as many clear pulsators as possible, rather than trying 
to push the limit on the pulsation amplitudes as low as pos- 
sible. Thus, most of our detected variables have rather high 
pulsation amplitudes, and the sample may still contain ob- 
jects with undetected amplitudes below 2 mma. Due to the 
shortness of our typical photometric sequences, there may 
also be multiperiodic pulsators among our null-detections, 
unidentified simply because their periods happened to be 
canceling each other out in a beat phase at the short time 
they were observed. The stars in Table [5] should therefore 
not be considered non-variable, rather, we prefer the term 
not observed to vary (NOV). 



For the aims of this campaign we consider a star as 
NOV if there is no peak in its Fourier Transform (FT) 
with an amplitude higher than 3.0 times the mean ampli- 
tude level (a). Howev er, we require an amplitude of 3.7 a 
(jKuschnig et al.|[l997l 99% confidence level) in order to call 
the object a pulsator. Objects that fall between these two 
tiers are considered pulsator candidates until further obser- 
vations can confirm the detection. Most pulsator candidates 
were reobserved in subsequent runs, but 11 objects remain 
in Table [5] with this classification, and will be discussed 
separately in Sect. 14.11 

Only the frequency region 1.5-10 mHz was used to com- 
pute a homogeneously. The low frequency range was ex- 
cluded to avoid the higher mean amplitude level caused by 
small trends in the light-curve due to extinction effects or 
sky transparency variations. 

3.1. PMT photometry 

The first observation run (Nl) was made with the Troms0- 
Texas 3-channel photoelectric photometer, equipped with 
Hamamatsu R647 photomultipliers. The data were reduced 
with the standard quilt softwar e package develope d for 
the wet (Whole Earth Telescope; leather et al1ll990l ). The 
reductions include smoothing of the sky background, and 
correction for extinction and linear trends. A total of 13 
stars were observed, and one pulsator detected. However, 
due to problems with the guiding many of the sequences 
were rather poor, giving a it of about 2 mma in the FT. 
Several stars showed a behaviour that could be interpreted 
as pulsations, but was actually due to periodic motion in 
the apertures. For this reason three stars were reobserved 
during the N2 run, and only the sequence with the lowest 
(7 is listed in Table [5l 

3.2. CCD photometry 

Further not observations were done with the High Reso- 
lution Adaptive Camera (hirac) for the second run (N2 
in Table [lj , and the Andalucia Faint Object Spectrograph 
and Camera (alfosc) for all the remaining runs. Both 
hirac and ALFOSC were equipped with Loral, Lesser 
thinned, 2048x2048 CCD chips, up until January 2004, 
when ALFOSC was upgraded with an E2V CCD 42-40 de- 
vice. The sky area available for locating reference stars 
was limited to ^3.7x3.7 and ^6.5x6.5 arcminutes 2 re- 
spectively for each of the two cameras. All CCD obser- 
vations used the multi-win dowed fast CCD photometry 
(|0stensen fc Solheiml [2000). Efficient fast photometry is 
easily achieved with this mode, since only small regions 
(windows) of the chip are read. 

Observation s were made with either a Bessell B-band 
filter (not #74. lBesseIilll990D . for the brighter targets, and 
a much wider filter (not #92) for the fainter ones 0This 
filter (hereafter referred to as the W, or wide filter) has 
the same centre as that of the V-filter (5 500 A) , but is a 
full 2 750 A wide, effectively covering all bands from B to 
R with more than 90% transmission. The cycle times were 
mostly set to 20 s, except for some of the brightest objects, 
which can have cycle times as short as 5 s, and the faintest 
ones, which may have times up to 40 s. Between 2.5 and 



4 Filter details with transmission curves for all NOT filters can 
be found on their web-pages; http://www.not.iac.es/ 
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Table 5. Truncated version of the catalogue of sdB stars observed during our search programme. 



Target Name 


Other Name 


RA 


Dec 


Wlntr 
PS 


Run 




At 


a 




Var 






J2000 


J2000 


mag 






s 


mma 


mma 




On 1-7 

bB 7 


MCI 0000-1637 


00:1)0:24.37 


— lb:21:0b.o 


12.7 


bo 


60 


20 


0.60 


1.26 


JNUV 


rlxL OOOz— zo4o 


Mtl 0002— 2o48 


00:Uo:0y.oo 


— zb:ol:4o.4 


15.6 


bo 


Q A 

39 


OA 

20 


1 A1 

t.ut 


2.U/ 


JNUV 


1 1 I ' Ann A 0707 

HJi 0004-2737 


TV/T^O^ AAA/1 0707 

Mtl 0004-2737 


00:06:46.26 


07 or» co /i 

—27:20:53.4 


13.7 


JN4 


60 


20 


0.44 


0.95 


JNUV 


ur nnn7 0010 
rirj U00 (—2212 


MO 1 000 1— 2 f 3 ( 


AA. fin. m 

U0:09:4o.91 


—21:56:14.4 


1 /I A 

14.9 


c ~ 

DO 


43 


1 n 
1U 


A TA 
0. (9 


1 AO 

1.48 


JNUV 


T\ T A 1 r I '' AAAO 1 A A C 

MCI 0008-1245 




AA11 OCA/1 

00:tt:25.04 


1 O OO A C A 

— 12:28:45.9 


14.3 


1 

b4 


93 


15 


0.59 


1.04 


JNUV 


KUV 04233+1502 




04:26:08.00 


+ 15:08:29.01 


14.3 


S4 


129 


20 


0.33 


1.01 




KPD 0716+0258 




07:18:57.85 


+02:53:14.67 


14.7 


T4 


270 


30 


0.25 


0.78 




KUV 07301+3659 




07:33:25.62 


+36:52:57.83 


14.8 


S4 


134 


18 


0.30 


0.93 




PG 1248+164 




12:50:50.23 


+16:10:03.40 


14.0 


T7 


65 


20 


5.07 


16.61 




PG 1313+132 




13:15:58.16 


+ 12:57:40.42 


14.1 


S6 


96 


20 


0.95 


3.35 




J1351+0234 




13:51:40.69 


+02:34:29.21 


17.1 


N7 


70 


30 


0.91 


2.93 




PG 1722+286 




17:24:11.97 


+28:35:26.90 


12.8 


Nl 


1191 


5 


1.01 


3.08 




PG 1725+252 




17:27:57.39 


+25:08:35.70 


12.7 


T2 


61 


10 


3.38 


12.39 




HS 2029+0301 




20:31:52.63 


+03:11:17.39 


15.6 


Nl 


952 


5 


2.13 


6.41 




KPD 2215+5037 




22:17:20.62 


+50:52:58.09 


13.6 


N5 


45 


25 


0.69 


2.39 





The full table with 285 entries is availa ble in electron ic form from CDS, via anonymous ftp to cdsarc.u-strasbg.fr or via 
http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A+A/???/??? . Only the first five entries in the table are listed, plus the ten stars 
that have peaks between 3cr and the 3.7 a cutoff required for the NOV flag (discussed in Sect. fO]) . The stars in the sample for 
which pulsations were detected are listed in Table[9] The columns include target names, coordinates, a magnitude estimate (usually 
photographic), the run ID from Table [U the number of data points in the light-curve (N p ) used to compute the FT, the sampling 
interval (At), the mean amplitude in the FT (a), the amplitude of the highest peak in the FT (j4 max ), and a non- variability flag. 



6.5 s of the cycle times are readout and shutter overheads, 
the exact number depending on the number of reference 
stars selected and the size of the readout windows. 

The CCD data were reduced on-line with the Real Time 
Photometry (RTP) program develop ed by one of us as part 
of his Ph.D-project (|0stensenll2OOQl see also Paper II). This 
software was particularly useful for the purposes of this sur- 
vey, as it performs real-time reduction of the data during 
sequencing, displays the light-curve of the target, and com- 
putes the FT. Thus, the noise level and the presence of clear 
variability could be checked at the telescope, although care- 
ful data reduction was also done afterwards. The processing 
includes bias level removal, flat fielding, sky subtraction, 
extinction correction and aperture photometry using aper- 
tures that track each star's geometrical centre. 

The optimal aperture for each sequence was selected 
after processing all data sets with apertures of a wide range 
of diameters and choosing the one that gave the lowest 
noise in the FT. The apertures tested ranged between 10 
and 40 pixels in radius, which corresponds to 1.9 and 7.6 
arcseconds on the sky, while the best results usually stayed 
within the range of 18 to 22 pixels (3.4 to 4.2 arcseconds), 
but on nights with particularly poor seeing the optimal 
radius could be as high as 30 pixels (5.7 arcseconds). 

4. Results 

A total of 309 stars deemed to lie in the predicted instabil- 
ity region for sdB stars were checked for variability during 
25 observing runs from 1999 to 2009. We did not apply 
a specific cutoff for T g and log g since our estimate of 
these parameters were preliminary and sometimes crude. 
Instead we prioritised our targets by their distance from 
the centre of the instability region at 34 kK. In the later 
runs we shifted our focus more towards the region where 
the hybrid pulsators are located. Of all our targets, 20 were 



observed to pulsate with frequencies corresponding to the 
V361 Hya stars, which corresponds to a ~ 6.5% ratio of 
fast pulsators/NOV sdBs. A further four stars were found 
to have such strong reflection effects that the slope was ob- 
vious even in our short light-curves. This fraction of 1.3% 
is only a lower limit for the sdB+dM population, as there 
may be quite a few more objects with longer periods, which 
our survey for short period pulsators was not sensitive to. 
In fact, we know of one more reflection variable in our sur- 
vey, BPS CS 22169-1, which we found to be constant in our 
half-hour light-curve. Since these results were not obtained 
at the NOT as part of this survey, they will be reported in a 
future paper. The 285 stars for which no obvious variability 
was detected are presented in Table [5] 

A different detection threshold was achieved for each 
target, depending on its magnitude, the filter used for the 
observations (B or W) , the length of the photometric time- 
series, and the weather conditions in each particular case. 
In the top panel of Fig. [T] the mean amplitude level in the 
FT (a) obtained for the sample is shown as function of each 
star's magnitude. A general trend of larger a for fainter ob- 
jects can be distinguished, although in most cases the time- 
series sequence was stopped when the real-time FT showed 
no significant peaks in the relevant region and the highest 
peaks reached amplitudes ~ 1 mma. Theoretical a levels as 
a function of magnitude for light-curves with 1, 5, 20, 100 
and 300 data points are also plotted in this panel. They 
represent the contribution to a of both the scintillation and 
photon shot noise, and thus should be understood as the 
limiting accuracy in the FT obtainable from a time series 
with the same number o f data points. The scintillation noise 
was computed following I Youngj (|1993h . and the photon shot 
noise was estimated using the signal-to-noise calculator of 
the not. In both typical airmass X = 1.2 and a 

typical exposure time of 20 s were used. Most points follow 
the predicted trends for light-curves with between 20 and 
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Fig. 1. Top panel: Photographic B magnitudes of the tar- 
gets as function of ct. The labeled curves indicate the the- 
oretical noise limits for a typical photometric observation 
with our system (1), and scaled to show the limit achieved 
for a time-series of 5, 20, 100 and 300 measurements. 
Middle: Maximum amplitude detected in the FT for each 
target, as function of ct. The pulsating sdBs from this sur- 
vey are marked with star symbols (published earlier) and 
diamonds (presented in this paper). The continuous line 
indicates the 3.7a level (99% confidence), and the dashed 
line indicates the 3. Oct level. Bottom: Histogram of the noise 
level achieved for the sample. 



300 points, and the outliers are mostly due to poor weather 
conditions. 

The middle panel of Fig. Q] shows, again as function of ct, 
the maximum amplitude obtained in the FT for every ob- 
served target within the 1.5-10 mHz frequency range. The 
continuous line indicates the 3.7ct (99%) threshold, mean- 
ing that if A max is above this line, the star is considered as 
pulsator. The dashed line indicates the 3. Oct (90%) thresh- 
old, and stars whose A max falls below this line are consid- 
ered as NOV, according to our definition. Fourteen of the 
twenty sdBVs discovered in this survey are marked with 
stars (published earlier) and four with diamonds (the new 
ones presented here); the remaining two pulsators, EPPsc 
and QQ Vir, have such large pulsation amplitudes that they 
lie beyond the range shown in the figure. Note that we have 
used only the discovery data for these plots, so the lowest 
amplitude pulsator, LM Dra, and the one on the 99% confi- 
dence line at ct = 0.75, PG 1657+416, have been confirmed 



in follow-up runs with higher confidence levels. The objects 
between the 3.0 and 3.7 ct confidence lines are discussed 
below, and the pulsators marked with diamonds will be 
discussed in Sect. [SI 

The bottom panel of Fig. [T] represents an histogram of 
the mean amplitude level achieved along the survey. For 
most of the targets, the noise level in the FT was below 
lmma (239 targets = 83%), and in 50% of the cases a 
ct < 0.6 mma was obtained. 



4.1. Pulsator candidates 

Ten objects in Table [5J have amplitude peaks between 3.0 
and 3.7ct, higher than our NOV requirement. They must 
therefore still be considered pulsator candidates. Three of 
these have a ct of 0.33 mma or less, and have been ignored. 
The others deserve particular mention. 

One object, PG 1248+164 has such a high ct that it is 
off th e scale in Fig. CD Th e object comes from the sam- 
ple of iMaxted et al.l ([2001) an d deserves further observa- 
tions. The object on the 99% confidence line at ct = 3.4 is 
PG 1725+252. Normally it would have been reobserved, but 
as it was real i sed th at this object is also in the survey of 
IBilleres et all (|2002ft with a comfortable 0.03% limit, no 
follow-up was performed. The object just below the 99% 
limit at ct = 0.7 is KPD 2215+5037, observ ed during run 
N5. T his star was also part of the survey of Bil leres et al.l 
(2002) with a limit of 0.2%. We note a rising trend in our 
light-curve on the 1% level over the 20 minutes covered by 
our observations, which could indicate that the object is a 
long period variable. 

The three objects between the 3ct and 3.7 ct confi- 
dence lines at -Ict are PG 1313+132, PG 1722+286 and 
J1351+0234. PG 1313+132 shows variations at low frequen- 
cies, consistent with g-mode pulsations, which is likely con- 
sidering its low temperature (T e g = 25.6kK, log<7=5.41; 
Saffe r et al1ll994L but uncertain as the observations were 
done at relatively high airmass. If we compute the noise 
level only for frequencies higher than 5 mHz we get ct = 0.53 
and the highest peak is well below three times this value. 
The light-curve of PG 1722+286 was obtained during the 
problematic Nl run. An attempt was made to reobserve it 
during T2, but poor weather resulted in an even worse light- 
curve. After a limit on pulsation of 0.08% was presented by 
IBilleres et al.l (|2002[ ). no further attempts were made to ob- 
serve this target. The last of the three, J1351+0234, is just 
noisy due to its low brightness (g' = 17.1). 

The three objects that lie on the 3ct confidence 
line around ct = 2 are HS 1813+7247, EGGR491, and 
HS 2029+0301 (only the last of these is actually above the 
3ct line, but we will discuss all three). HS 1813+7247 has a 
slope in the light-curve, but it is not consistent between the 
different reference stars, and so is most likely due to differ- 
ential extinction as the W-band filter was used. EGGR491 
is another faint target from the SDSS sample, so the light- 
curve is just noisy, with no obvious features. HS 2029+0301 
is also from the Nl photoelectric run, and has not been 
reobserved. Further observations would be required to give 
a more useful limit on any pulsational behaviour of this 
object. 
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Fig. 2. Amplitude spectrum of HE 2151-1001. The lower 
panel shows the residual spectrum after prewhitening the 
two frequencies listed in Table [7] from the light-curve. The 
dashed line indicates 3.7 times the noise level. 

5. New pulsators 

The four pulsators found in our sample and not published 
elsewhere are HE 2151-1001, HS 2125+1105, PG 1033+201, 
and HE 1450-0957. The details of the discovery observa- 
tions as well as follow-up observations on the first two done 
with the NOT are provided here. Table [5] includes the de- 
tailed log of observations for all these new sdBVs. 

For all the pulsators we derived the frequency content 
of the light-curves by performing a non-linear least-squares 
fit to a sine function, using the highest amplitude frequency 
in the amplitude spectrum (top panels of Fig. [2j [3j |4] and 
[3]). This fit was then subtracted from the light-curve to 
compute the residual amplitude spectrum (bottom pan- 
els). Notice that the noise level is recomputed after each 
prewhitening step, so that the 3.7cr level lines in the figures 
drop after removing the established signal. If there was still 
any peak left above the significance threshold, the original 
light-curve was fitted to the sum of two sine functions with 
the established frequencies. Table[7]lists the best fitting fre- 
quencies and amplitudes, with the associated errors from 
the least-squares procedure. 

5.1. HE 2151-1001 

This target was observed for ^40min during the run S2. 
Its amplitude spectrum is shown in Fig. [51 where some sig- 
nal above the threshold can be distinguished at ^8 mHz. 
This frequency range is typical for a pulsating sdB with 
these T c ff-log<? parameters (sec Tabled]). 

Our non-linear least-squares frequency analysis revealed 
two low amplitude (3mma) peaks, as listed in Table [7] 
Note that the frequency distance between the two peaks 
(~140/HIz) is within the frequency resolution (^325 /iHz). 
Another short sequence on HE 2151-1001 was obtained dur- 
ing S3 (see Table , but the noise level achieved was not 
sufficient to detect the oscillations above the significance 
threshold. 

This target was observed spectroscopically as a white 
dwarf ca ndidate by the SPY survey. The spectrum was anal- 
ysed by iLisker et al.l (|2005l ). who found no spectral indi- 
cation of a companion, but noted that it has a peculiar 
Ha profile. The 2MASS photometry (J =16.53, H= 16.37, 
K= 15.91) indicates a main sequence companion, but 
this is uncertain as the values are close to the faint 
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Fig. 3. Amplitude spectrum of HS 2125+1105. The lower 
panel shows residual spectrum after subtracting one fre- 
quency from the original light-curve. The dashed line in- 
dicates 3.7 times the noise level. The peak at 7.3 mHz is 
just above the 3.7 a level if the low frequency domain is 
excluded when estimating a. 



limit of the 2MASS survey. The NOMAD c atalogue (Naval 
Obse rvatory Merged Astrometric Dataset. IZacharias et al.l 
l2005h provides reasonably good photographic magnitudes 
(5 = 15.25, V = 15.57, R = 16.09), and the UV magnitudes 
FUV = 14.55 NU V =14.87) from the galex satellite 
Martin et al.l I2005T) all indicate colours typical for a hot 
subdwarf. 

5.2. HS2125+1105 

The faint (B = 16.3) HS 2125+1105 was first observed dur- 
ing run S2. A 6.4 mma peak was detected at 6.89 mHz, 
although its significance level was exactly 3.7 times the 
noise level (see Table [5]). A year later, a slightly longer 
light-curve with better photometric conditions was ob- 
tained during run ST, which reduced the noise level to 
0.95 mma. However, a peak detected at the same frequency 
also dropped its amplitude to 3.6 mma, again at exactly the 



Table 6. Details from the discovery runs for the four new 
sdBVs. 



Target Run N p At a A max 

ID s mma mma 



HE 2151-1001 


S2 


120 


30 


0.58 


2.20 




S3 


73 


30 


1.05 


2.35 


HS 2125+1105 


S2 


67 


35 


1.75 


6.45 




ST 


95 


30 


0.95 


3.56 




S3 


116 


42 


0.72 


3.42 




S3 


433 


30 


0.40 


4.35 




S3 


541 


30 


0.33 


3.84 




S4 


89 


28 


1.17 


4.76 


PG 1033+201 


T7 


118 


20 


1.01 


4.95 


HE 1450-0957 


S6 


125 


30 


0.90 


6.81 



'Run ID' refers to the identifiers listed in Table [TJ N p is the 
number of points for the specific sequence, At is the sampling 
time interval, a is the mean noise level in the FT, and A max is 
the amplitude of the highest peak. 
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Fig. 4. Amplitude spectrum of PG 1033+201. The lower 
panel shows the residual spectrum after prewhitening the 
two frequencies listed in Table [7] from the light-curve. We 
have also plotted the amplitude spectrum of the first half 
(red dashed curve) and second half (blue short dash, long 
dash curve), to demonstrate that the main peak is present 
with the same amplitude in both halves of the dataset. The 
dashed line indicates 3.7 times the noise level. 



3.7 times the noise level. During run S3, the object was ob- 
served on each of the three nights as indicated in Table |H1 
and the good conditions finally allowed us to get reliable 
confirmation of the main peak. The amplitude spectrum 
obtained from the combined light-curve is shown in the up- 
per panel of Fig. [21 In Table [7] the results of a frequency 
analysis are listed. A 4 mma peak at 6.83 mHz is detected, 
but this time at a reassuringly high S/N. In this sequence 
a second peak at 7.34 mHz was also detected, as listed in 
Table [7] The bottom panel of Fig. [3] shows the residual 
amplitude spectrum after prewhitening the main frequency 
from the light-curve. 

In order to check the stability of this mode, an addi- 
tional short sequence was taken in 54. The same /i was 
found, this time with a 4.7 mma amplitude. 

This H S star was not par t of the original sample of sdB 
stars from lEdelmannl (|2003lh but was observed as part o f 
the SPY survey and was analysed by iLisker et al.l |2005) . 
The star has reliable optical photometry (but no spec- 
trum) in the SDSS survey (u' = 16.13, gf = 16.36, r' = 16.78, 
i' = 17.09, z' — 17.33). There is no indication of a main se- 
quence companion from the spectroscopy or SDSS photom- 
etry. The object is well below the faint limit of the 2MASS 
survey, so no IR colours are available. The GALEX photom- 
etry is FUV=15.6, NUV = 15.8. 

5.3. PG 1033+201 

PG 1033+201 was one among 15 targets observed with the 
£>-band filter during the technical night T7. During the run 
we had some problems with the observing system, which 
resulted in some frames in the photometric sequence be- 
ing dropped. For this reason the real-time processing pro- 
gram was unable to process the entire sequence. Due to 
these problems we stayed on target almost 40 minutes, 
which is considerably longer than our standard sequence 
on such bright targets. Eventually, we terminated the ob- 
servations in order to reinitialise the CCD system, but we 
never returned to reobserve this star. During the careful 
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Fig. 5. Amplitude spectrum of HE 1450-0957 (top) and 
residual amplitude spectrum after subtracting the main 
peak (bottom), with 3.7ct noise levels (dashed lines). 



reprocessing of the photometry for this paper, we excluded 
two lost and one bad data point from the 118 frame se- 
quence, and performed a regular frequency analysis of the 
cleaned light-curve. The amplitude spectrum is shown in 
the top panel of Fig. 0] as a continuous line. A clear peak at 
5.8 mHz is well above the 3.7<r (horizontal dashed) line. We 
divided the light-curve into two equal halves, and we can 
see from the corresponding spectra (also shown in Fig. H] 
as red dotted and blue dash-dotted curves), that the main 
period is present in both halves with the same amplitude. 
A second frequency at 6.8 mHz is significant only in the 
first part, indicating that it is either spurious or a com- 
bination of frequencies suffering cancellation from beat- 
ing within our short sequences. This exercise demonstrates 
how short-period pulsations can be detected with reason- 
able confidence even in short (^20 m) light-curves, as long 
as one acquires adequate signal-to-noise in the individual 
data-points. The bottom panel displays the residual ampli- 
tude spectrum after subtracting the frequency at 5.8 mHz. 
The remaining peak at ~6.8 mHz becomes significant when 
the noise is computed in the prewhitened spectrum (see 
Table 0). We conclude that PG 1033+201 appears to be 
a typical short-period sdB pulsator with a main pulsation 



Table 7. Pulsation properties for the four new pulsators. 



Name 


Id 


Frequency 


Amplitude 


S/N 






DiHz] 


[mma] 




HE 2151-1001 


h 


7931 ± 63 


3.04 + 1.81 


6.5 




h 


7791 ± 64 


3.00 + 1.80 


6.5 


HS 2125+1105 


h 


6835.66 + 0.07 


4.05 + 0.19 


19.6 




h 


7342.35 + 0.35 


0.85 + 0.19 


4.1 


PG 1033+201 


h 


5850 ± 27 


4.9 + 0.6 


8.1 




h 


6837 + 53 


2.5 + 0.6 


4.5 


HE 1450-0957 


h 


7202 ± 8 


6.9 + 0.4 


13.7 




h 


8475 ± 21 


2.5 + 0.4 


5.2 



Pulsation frequencies and amplitudes are given as derived by 
non- linear least-square fitting to the light-curves. The S/N is 
computed after prewhitening the detected frequencies from the 
signal. 
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amplitude of 5 mma, with likely further frequencies at the 
2 mma level. 

This object w as included in our sample based on the 
decomposition by lAllard et al.l (|l994f) , who estimated the 
system to be an sdB+F9 composite from their Cousins 
BVRI photometry (V = 15.67, B-V = -0.17, V-R = -0M, 
R-I = +0.06), and provided a colour temperature for the 
primary of 31.5 kK. The object occurs in the SDSS survey 
as J103638.93+195202.2, and has photometry (u 1 = 15.14, 
g' = 15.39, r' = 16.80, %' = 16.01, z' = 16.21) but no spec- 
troscopy, which prevents us from acquiring a better tem- 
perature and a gravity estimate. The 2MASS IR photome- 
try clearly supports the notion of a main-sequence compan- 
ion (J= 15.46, =15 16, K = 15.13, i.e. J-H> +0.3, see 
iReed fc Stienrnll2?5ol . 

5.4. HE 1450-0957 

HE 1450-0957 was one of the three targets surveyed during 
the final run of the programme (S6). A rather long light- 
curve was obtained (lh5min), since its oscillations were al- 
ready clear as the data were processed as they arrived. Its 
amplitude spectrum (upper panel of Fig. [5]) reveals the mul- 
timode behaviour of this target. At almost 7 mma, the main 
peak at 7.2 mHz has the highest amplitude of the four new 
pulsators presented here. A significant second peak is also 
detected in the amplitude spectrum of HE 1450-0957 (see 
bottom panel of Fig. 5), but more low amplitude modes 
are likely to remain in the range 6 -9 mHz. The frequen- 
cies derived from this discovery light-curve are also listed 
in Table 

This object c ame into our samp le from the HE stars 
surveyed by SPY (jLisker et al.l 120051 ). It is included in the 
Edinburgh-Cape survey as EC 14507-0957 where reliable 
photometry can be fou nd (V = 15.27, U-B = -1.04, B- 
V = -0.21). As noted bv lReed fc Stienind (|2004l ). its 2MASS 
magnitudes (J— 15.58, if = 15.41, K= 15.36) are indica- 
tive of a main sequence companio n, but this compan ion 
must be very weak in the optical as [Lisker et al.l (|2005f) do 
not detect any trace of it in their VLT/uves spectrum. 



6. Summary and discussion 

The NOT search programme for pulsating sdB stars has 
been a great success, contributing twenty new pulsators to 
the known population of short perio d sdBVs. Together wit h 



(2009) 



2010) 



trie two new pulsators presente d byUv nkcnny et al. ( 
another two recently found by iBarlow et al.l ( 200E , _ 
and HE 0230-43231 our latest tally brings the total up to 
49. Our effort to implement an efficient system for CCD 
photometry at the NOT has certainly been fruitful, produc- 
ing on average one new pulsator for every three nights of 
observations. This is a substantially higher efficiency than 
any other group has reported. 

6.1. Follow-up 

Many of our pulsators have already received detailed follow- 
up. Most effort has been put in on V391 Peg, which was the 
subject of several campaigns, and its pulsation periods are 



still being monitored regularly (|Silvotti et al.l |2007[ ) . The 
original aim of this monitoring was to detect evolutionary 
period changes, but the periodic variation that was found 
is more likely to be caused by an orbiting planet induc- 
ing period variations through the light travel time effect. 
Several other stars in our sample are now regularly moni- 
tored for vari ations in the puls ation period under the exo- 
time project (|Lutz et a l. 2009). As a substantial time base 
of these observations develops, we expect to detect more 
planet-hosting hot subdwarf stars. 

QQVir was the first sdBV for which spectroscopic 
line-profile varia tions were used to cons train its main 
pulsation mode ()Telting fc Ostensenl |2004|) , and was the 
subject of an exte nsive photometric campaign in 2003 
(|Silvotti et al.|[2006l) revealing 15 pulsation frequencies. The 
detailed frequency spe ctrum from this campaign allowed 
ICharpinet et al.l ((2006) to obtain an asteroseismic solution 
using the forward method. QQVir also became the first 
star from our sample to be targete d with high-resolut ion 



spectroscopy with the VLT in 2008 ( Telting et al 



LMDra was followed extensively bv lReed et al 



2010ft. 



(I2007aft 

with multisite observations between 2003 and 2004, and in 
spite of its extremely low pulsation amplitudes they de- 
tected six pulsation periods with amplitudes between 1.0 
and 2.2 mma. V42 9And and V1636 0ri were followed by 
IReed et all (|2007bt) in Nov-Dec2005, revealing a rich fre- 
quency spectrum with fourteen significant peaks between 
0.45 and 4.35 mma in V429And, and a simpler spectrum 
with only thre e significant peaks b etween 0.8 and 11.1 mma 
in V1636 Ori. IReed et al.l (|2007bh also report that LS Dra 
appears to be monopcriodic, but with a highly variable am- 
plitude, changing from a maximum of 5.26 to a minimum of 
0.87 mma over their campaign spanning 47 days. They fur- 
ther confirm three of the four pulsation periods in V387 Peg. 

J1717+5805 was targeted with high-speed simultaneous 
3-channel p hotometry u s ing U LTRACAM on the WHT, as 
reported by 
described in 



Aerts et a l. ( 2006), resolving the broad peak 



Paper vl into two separate frequencies. 



5 Originally r eported as "an unusual hot subdwarf pulsator" 
by iKoenl {2007) but now confirmed to be a more regular sdBV 
by Kilkenny (priv.comm) 



6.2. Summary 

The four new pulsators presented in this paper certainly 
deserve more attention as well. Confirming the pulsations 
in PC 1033+201 would be a priority, but an intensive cam- 
paign would probably be required to disentangle its multi- 
periodic nature. The pulsations observed in HE 2151-1001 
are barely above our significance threshold and need to be 
confirmed. HS 2125+1105 is clearly a stable pulsator, as the 
same oscillation period was detected in six sequences span- 
ning four years, but we would suspect more frequencies to 
be discovered with more dedicated follow-up. HE 1450-0957 
is also a convincing case where residual power in the ampli- 
tude spectrum between 6 and 10 mHz indicates the presence 
of unresolved pulsations that could be revealed with ded- 
icated follow-up. The existence of pulsations in this star 
was recently confirmed by observations from South Africa 
(C. Koen, priv. comm.). 

In Table[S]we compare the fraction of pulsators detected 
in each of the different samples we have drawn our candi- 
dates from. A fraction of pulsators between 5 and 12% is 
clear when targets are selected based on spectroscopic tem- 
perature estimates. We estimate that it should be possible 
to obtain a pulsation fraction of about 10 %, when starting 
with an undepleted sample and reaching a noise level below 
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Temperature [kK] Period [s] 

Fig. 6. a) T e ff-logg diagram for the objects in Tabled Filled circles indicate V361Hya type pulsators and diamonds 
indicate hybrid DWLyn type pulsators. The size of the symbols is proportional to the pulsation amplitude, A max . The 
nineteen sdB Vs from this survey that we have T e s and log g for are shown in blue, and short period pulsators from other 
surveys in red. The green lines indi cate the canonical zero age EHB (lower line) and terminal age EHB for an 0.47 M© 
core ffrom lKawaler fc Hostlcr 2005) with envelope mass increasing from bottom to top. The ZAEHB starts at a logg of 
5.9 with an envelope mass fraction of only 0.02%, and the objects that appear to lie on the extension of this line cannot 
be explained by the canonical models, b) P-log g diagram for the same pulsators. The bars indicate the range of detected 
pulsation periods in any particular star, excluding harmonics and g- modes. The log g values have been shifted by up to 
±0.04 dex to avoid overlaps. 



1 mma. This could probably be increased by focusing only 
on the hot end of the EHB, between 31 and 36 kK, at the 
cost of sacrificing discoveries of the rarer pulsators at the 
cool end of the instability region. 

That the colour selected sample has a lower success frac- 
tion than the samples based on spectroscopy is not sur- 
prising, but it is su rprising that there are so few pulsators 
in the sample from ISaffer et al.l ()1994l ) . The total number 
of sdB stars in that sample hotter than 28 kK are 44, of 
which we observed 16 (plus another eight that were out- 
side the borders of the instability region). Excluding the 
known g-mode pulsators on the cool end of that sample, 
we have two short period sdBVs in 41 stars; QQ Vir fo und 
by us and V1078Her found by iBonanno et all (|2003T >. A 
large number of sta rs from this sample were surveyed by 
iBilleres et all (|2002h . and the remaining stars have been 
checked by other observers (but no limits have been pub- 
lished), so we do not expect more clear pulsators in this 
sample, although some low amplitude variables may have 
been missed. Similarly, the n umber of sdB s t ars in the insta- 
bility region in the sample of lMaxted et al.l (|200ll ) is 34. We 
found no pulsators in this sample, but it has been heavily 
exploited by other groups; [Billeres et al.l (|2002D surveyed 14 
of them and found one: UYSex. When combining informa- 
tion from this and other unpublished surveys, we can only 
find a single star tha t has yet to b e surv eyed. Thus, the 
best fraction from the iMaxted et al.l (|200ll ) sample we can 



infer is 3 %, which is lower than what we have found in the 
HS, SPY and SDSS samples, but with such small samples 
this is barely significant. 

6.3. The big picture 

We have compiled all the known short period sdB pulsators 
in Table IH1 using data from the literature and added the re- 
sults from this paper. We included prelimin ary physical pa- 
ramet ers on two new pulsators published bv [Kilkenny et all 
(2009), based on recent low-resolution spectroscopy from 
not/alfosc (3250-6150A at — 4.5A resolution), using 
the same atmosphere models as in our other determina- 
tions. A few stars still lack temperature and gravity esti- 
mates from spectroscopy, but we have supplied values from 
our own ALFOSC spectroscopy where available. The table 
may also not be entir ely complete. PG 0856+121, for which 
iPiccioni et al.l ()2000D reported what they termed 'a pul- 
sation episode', was dropped from the table as repeated 
attempts (Reed, priv. comm) have failed to confirm any 
pulsations. Such brief pulsation episodes may be common 
in sdB stars, but they are hard to catch and characterise. 

The first sdB p ulsator candidate in a globular clus- 
ter was reported bv iRandall et al.l (|2009al ). but while the 
single clear period at 114 s makes this star a very likely 
sdBV, spectroscopic confirmation remains to be done. The 
same is the case for the pulsator candidate reported by 
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Table 9. All short period sdBV stars published to date. 



Name 


Survey name 


mv 


T e ff 


log 9 


logy 


C 


Pi 


angc 


A m ax 


N P 


References 






mag 


kK 


dex 


dex 






S 


mma 






V361 Hya 


EC 14026-2647 


15.3 


34.7 


6.10 




G2 


134- 


-144 


12 


2 


1,4 


EO Cet 


PB 8783 


12.3 


35.7 


5.70 




F0 


94- 


-136 


9 


7 


2.4 


UXSex 


EC 10228-0905 


15.9 


33.5 


6.00 




G2 


140- 


-152 


14 


3 


3.4 


V4640 Sgr 


EC 20117-4014 


12.5 


34.8 


5.87 




F5 


137- 


-158 


4 


3 


4 


UYSex 


PG 1047+003 


13.5 


33.2 


5.80 


-2.00 


ah 


104- 


-175 


10 


18 


5,fel 


NYVir 


PG 1336-018 


13.5 


31.3 


5.60 


-2.93 


M5 


97- 


-205 


10 


28 


6,fe2,u2 


V2203 Cyg 


KPD 2109+4401 


13.4 


31.8 


5.79 


-2.23 


a 


182- 


-213 


9 


8 


7,zl,hl 


V338 Ser 


PG 1605+072 


12.8 


32.3 


5.25 


-2.53 


a 


350- 


-573 


64 


50 


8,12 


KL UMa 


Fcigc 48 


13.5 


29.5 


5.50 


-2.93 


WD 


343- 


-379 


6 


4 


9,hl,o2 


DTLyn 


PG 0911+456 


14.6 


31.9 


5.77 


-2.55 


a 


149- 


-192 


7 


7 


10, r5 


KY UMa 


PG 1219+534 


13.2 


34.3 


5.95 


-1.50 


a b 


128- 


-149 


g 


4 


10, hi 


AQ Col 


EC 05217-3914 


15.6 


31.3 


5.76 


-3.00 




215- 


-218 


8 


8 


H,r3 


V1405 Ori 


KUV 04421+1416 


15.1 


32.0 


5.72 


-2.50 


dM 


185- 


-232 


20 


7 


11,7-4 


V2214 Cyg 


KPD 1930+2752 


13.8 


35.2 


5.61 


-1.50 


WD 


145- 


-332 


5 


44 


13, gl 


LMDra 


PG 1618+563B 


13.5 


33.9 


5.80 


-1.60 


F3 C 


108- 


-144 


6 


2 


14, r3 


DV Lyn 


HS 0815+4243 


16.1 


33.7 


5.95 


-2.10 




126- 


-131 


7 


2 


15 


V384 Peg 


HS 2149+0847 


16.5 


35.6 


5.90 


-1.80 




142- 


-159 


11 


5 


15 


V391 Peg 


HS 2201+2610 


13.6 


29.3 


5.40 


-3.00 


PI 


347- 


-351 


10 


3 


15, si 


LSDra 


HS 1824+5745 


15.6 


33.1 


6.00 


-1.52 






139 


5 


1 


16, rl 


V387 Peg 


HS 2151+0857 


16.5 


34.5 


6.10 


-1.37 


F-K 


129- 


-151 


4 


5 


16, rl 


V429 And 


HS 0039+4302 


15.1 


32.4 


5.70 


-2.20 


a 


182- 


-235 


8 


4 


16 


V1636 0H 


HS 0444+0458 


15.4 


33.8 


5.60 


-1.85 


a 


137- 


-168 


13 


2 


16 


EKPsc 


PG 0014+067 


16.3 


34.1 


5.77 


-1.68 


b 


137- 


-169 


3 


13 


17,cl,ul 


EPPsc 


PG 2303+019 


16.2 


35.3 


5.74 


-1.70 


MS? 


128- 


-145 


16 


3 


18 


QQVir 


PG 1325+101 


13.8 


34.8 


5.81 


-1.65 


o. b 


94- 


-168 


26 


14 


18,tl,s2 


DWLyn 


HS 0702+6043 


14.7 


28.4 


5.35 


-2.70 


a 


363- 


-384 


22 


2 


19, s3 


V1078Her 


PG 1613+426 


14.1 


34.4 


5.97 


-1.65 


a 




144 


5 


1 


20 


PG 0048+091 




14.3 


34.2 


5.69 


-3.00 


F5 


90- 


-192 


6 


28 


21,r3,t 


PC 0154+182 




15.3 


35.8 


5.80 


-1.67 


G-K 


110- 


-164 


10 


6 


21,7-1,} 


J1717+5805 




16.7 


34.4 


5.75 


-1.80 


G9 


137- 


-144 


6 


2 


22, f 


Balloon 090100001 




11.8 


29.4 


5.33 


-2.54 


ab 


178- 


-356 


60 


3 


23, ol 


PG 1419+081 




14.9 


33.3 


5.85 


-1.80 


G5 




143 


7 


1 


24, t 


J1445+0002 




17.4 


35.9 


5.75 


-1.60 


GO 


120- 


-142 


8 


3 


24, t 


J1642+4252 




15.7 


32.4 


5.80 


-2.00 


GO 


130- 


-138 


3 


2 


24, t 


EC 09582-1137 


PG 0958-116 


15.2 


34.8 


5.79 


-1.68 


a 


136- 


-169 


8 


5 


25, r6 


EC 11583-2708 




14.4 








MS 


114- 


-149 


3 


4 


25 


EC 20338-1925 


BPS CS 22880-18 


15.6 


35.5 


5.75 


-1.71 


a 


135- 


-168 


27 


5 


25,} 


RATJ0455+1305 




17.2 


29.2 


5.20 






363- 


-374 


19 


2 


26,61 


PG 1657+416 




15.9 


32.2 


5.80 


-1.60 


G5 


125- 


-143 


4 


5 


27,} 


2M0415+0154 




14.0 


34.0 


5.80 


-1.60 


a b 


144- 


-149 


6 


3 


28 


EC 01541-1409 




12.2 


37.1 


5.71 


-1.21 


a 


78- 


-146 


10 


6 


29,} 


EC 22221-3152 




13.4 


35.1 


5.85 


-1.47 


a 


84 


-176 


9 


10 


29,} 


JL166 




15.0 


34.4 


5.75 


-0.80 


MS 


97 


-167 


4 


10 


30 


CS 1246 


CS 124636.2-631549 


14.6 


28.5 


5.46 


-2.00 


b 




372 


35 


1 


31 


HE 0230-4323 




13.8 


31.6 


5.60 


-2.58 


dM 


282- 


-310 


8 


3 


32, 11 


HE 2151-1001 




15.6 


35.0 


5.70 


-1.60 


MS 


126- 


-128 


3 


2 


t,n 


HS 2125+1105 




16.3 


32.5 


5.76 


-1.86 


a 


136- 


-146 


4 


1 


t,!l 


PG 1033+201 




15.4 


31.5 






F9 


146- 


-171 


5 


2 


t,al 


HE 1450-0957 


EC 14507-0957 


15.3 


34.6 


5.79 


-1.29 


MS 


118- 


-139 


7 


3 


t-'l 



Pulsators discovered by our survey are listed with their names in bold face, and hybrid DW Lyn type pulsators in italics. C gives 
the class of the companion, where available, and the notes in that column mark: a No IR excess from 2MASS or SDSS photometry, 
6 No detectable RV variations on a time-scale of hours or days. c The F3 companion to LMDra is separated by 3.7", but it is 
unclear whether the subdwarf itself is single. P ra n g e gives the range of periods observed, excluding marginal detections, harmonics 
and g-modes, A max gives the highest amplitude reported, and Np the number of independent periods reported. The keys to the 
references in the last column are given on the next page. 



ISilvotti et al.l (|2009|) , which was found in the Kepler satel- We note that the number of pulsators in Tabic with 
lite field-of-view, with a likely pulsation period of 125 s F-K companions are 18 of 49 = 37%, which is compatible 
found at the 3.5cr level. with the fraction of spectroscopic composites we found in 



0stensen, Oreiro, Solheim et al.: A Survey for Pulsating sdB Stars 



15 



Table 8. Comparison with the results of other surveys for 
pulsating sdBs (upper part) and our own survey divided 
into subsamples (lower part). 



Sample 


Size 


Pulsators 


Fract. 


Preselection 




# 


S 


c 


% 


on T eB 


South Africa 


~1200 


13 


7 


1.7 


No 


Billeres 


74 





4 


5.4 


Yes 


Dreizler 


12 


1 





8 


Yes 


HS 


80 


9 





10 


Yes, -comp. 


HE 


23 











Yes 


SPY 


27 


3 





11 


Yes 


SDSS 


40 





5 


12.5 


Yes 


BG 


11 


1 





9 


Yes 


Saffer 


24 


1 





4 


Yes 


Maxted 


24 











Yes 


colour 


~50 





2 


4 


Poor 



'Size' indicates the number of targets checked for variability in 
each sample. The number of pulsators discovered are separated 
into apparently single (S) and spectroscopic composites (C), but 
note that binaries with invisible (WD or dM) companions are 
included in the column for single pulsators. We also specify if 
the sample was preselected to contain only subdwarfs in the 
instability region or not. Note that a few of the pulsators are 
present in more than one sample, and that we keep the HS and 
HE stars observed only as part of the SPY survey apart from the 
HS and HE samples. 



the SDSS as a whole. Pulsators also occur in systems with 
close white dwarf or M-dwarf companions, which means 
that pulsations appear independently of whether the subd- 
warf was formed through stable Roche-lobe overflow or via 
common envelope ejection. The stability of the pulsations 
for stars such as V391 Peg has revealed that some sdBs are 
definitely single, which implies that a formation channel 
must exist that produces single sdB stars that are other- 
wise indistinguishable from those produced via the channels 
that produces wide sdB+F-K or short period binaries. 

In Fig. [5^ we show the location in the (T c g , log g) plane 
of the pulsators from our sample (blue) together with pul- 
sators from the literature (red) , as listed in Table [9j with 
symbol size proportional to A max . As our sample stars are 
fairly evenly spaced along the canonical EHB strip, it is 
immediately clear that pulsations are more common on the 
hot side than on the cool side, but the amplitudes are often 
much higher in the cool, low gravity region. The figure does 
not really reproduce the gap around T c g ~ 30 kK between 
the main group of pulsators on the hot end of the EHB and 
th e cool group of q- mode and hybrid pulsators seen in Fig. 3 
of I0stens en (2009), where the sp ectroscopy was ob tained 
from the independent BG survey ([Green et al.ll2008h . This 
may be due to the fact that diverse methods and model 
grids are used in the temperature determinations collected 
from the literature, and this can produce scatter as high as 
2 kK in T e ff and 0.2 dex in log g. If we plot only the stars for 
which we ourselves have obtained the fits on a consistent 
model grid, the gap is indeed present. 



Refere nces to 

(|1997t). 2 Koen et al 
^lO'Donoghue et all 
O'Donoghue et aT| 
'j Billeres et all (119981 1 



Table [5] Discovery papers: ^Kilkenny et al.l 



3 IStobie et al 
5 Billeres et al 



Koen et alJ (ff998b). 



;i997 
(11997 

H998TI. ^Kilkenny et all 
iKoenl fi"99Sf} . B lKoen eTaLI (11998al), 
lu lKoen et ail (|1999bl ). ^ iKoen et alJ 



The main cluster of pulsators correlat es very well with 
the in itial instability strip predicted by ICharpinet et all 
(|1997l ). with only a few rare pulsators appearing to be 
anomalous. While the sample presented here contains no 
stars hotter than 39 kK, we have observed quite a few such 
stars without finding any with significant variability. But 
these are classified as sdO or He-sdO stars, and will be pre- 
sented in a future paper. Candidates in that sample were 
selected from very diverse surveys and ended up having a 
wide range of physical parameters, as stars matching the 
properties of the on ly known sdO pulsator, J16007+0748 
(jWoudt et aHl2006P are extremely rare, preventing an ef- 
fective preselection such as we made for the sdB sample. 

Figured shows the surface gravities plotted against the 
range of pulsation periods listed in Tabled As the logg val- 
ues are derived with many different methods, the system- 
atic errors can be substantially larger than ±0.1 dex. Still, 
the expected relationship between pulsation periods and 
surfac e gravity seems to hold up quite well fsee lKoen et al.l 
Il999bl ). The relationship appears to break down for objects 
with surface gravities higher than 5.8 dex, but this may not 
be a real effect. Evolutionary tracks for EHB stars with 
any envelope thick enough to sustain pulsations can hardly 
reach higher surface gravities than about 5.9 dex (note the 
lower end of the canonical ZAEHB in Fig. [5]). On the other 
hand, it is well known that the very high rotation rates in 
binary systems, high metallicity, as well as the pulsations 
themselves, give rise to broadening of atmospheric lines. If 
these broadening effects are not accounted for in the models 
used for estimating the atmospheric parameters, the fitting 
procedure may converge towards a too high gravity to com- 
pensate. 



7. Conclusions and outlook 



The 
the 



interest in sdB pulsators can 
future, when satellite missions 



only increase in 
such as Kepler 
(jChristensen-Dalseaard et ail 120071 ) start to find numer- 
ous new sdBVs and observe them with unprecedented fre- 
quency resolution and at lower pulsation amplitudes than 
have ever been achieved from the ground. As we enter 
the age of space based asteroseismology for the faint sub- 
dwarf stars, we can expect our understanding of the in- 
cidence of pulsations and their amplitude variability over 



(1999a), 12 
?4 | Paper J ! 
"iPaper I 
(2003), 



Kilkenny et~aTl (119991) 



iPaper ii 
3, 19 lDreizler et all 
21 lKoen et ail iEw 



13 lBilleres et all ([200 
Paper iij. 17 | Brassard et all d200 



(2002) 



20 



Bonanno et al 
lAerts et al 



Paper yj- 
Pap er vj. ^Kilkenny et al. 

7 Paper vil 2s lPaper vnA 2M 

30| Barlow et al.1 (|2009l ). 3l lBarlow etaP 



(20 06), 2 3 lOreiro et al.1 (200 
( 20 0). 26 lRamsav et al.1 1200 
Kilkenny et al.1 (|2009l) 

(|2uT0l ). 32 lKoenl (|2007l ). Kilkenny (priv.comm), 1 This 
paper. Supplementary data: al [Allard et al. ( 1994) 
bl lBaran fc Fox Machadol IpOlOl), cl |Charpinet eTaTl ([2 005) 
9l lGeier et all (pool). fe THeber et al.1 IpOOOl), fcl lKilkennv et al 
(I2002f). ^Kilkenny et all (120031), n lLisker et al.1 (120051) 



^lOreiro et al.l (120051) . ^ lO'Toole etHI (12004D 



HReed et al 
(20093) 



(2007a), 



^Reed et al 
^iReed et al l (l2007bl) 

rblr 



2007) 
20071 

" iTelting fc 0stensen 



' Randall et al 
s2 Silvotti et al 



Reed et al.l 
(120101), ^iRandall et al 
Silvotti et atl (l2002a| . 
(2006), SJ lSchuh et al.l (l20od ), 

hooi), ^IVuckovic et all (|2006l ). " 2 IVuckovic et al.1 (|2007l) 
zl lZhou et all (|2006l ). * Our determination of T eff , log g and 
logy. 
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time and across the instability re gion to improve, as MOST 
(| Walker et al.ll2003h and CoRoT ((Michel et al.ll2008h have 
done for brighter pulsating stars. Meanwhile, a new gener- 
ation of ground based instruments are coming on-line, ca- 
pable of high temporal resolution and simultaneous multi- 
colour photometry, by splitting the light with dichroics and 
reading out several CCDs simultaneously. ULTRACAM has 
pioneered this technology, permitting reliable amplitude ra- 
tios to be measured and used directly for mode identifica- 
tion. 

As theoretical models of the interior structure of these 
stars become sophisticat ed enough to accura tely predict 
the pulsation frequencies (|Hu et al.ll200ll2009h . new high- 
precision photometric observations will help us distinguish 
between the possible formation scenarios. When this hap- 
pens we will finally be in position to unravel the complete 
evolutionary history of the hot subdwarf stars. 
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